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Applying a macroscopic reduction procedure on the improved quantum molecular dynamics 
(ImQMD) model, the energy dependences of the nucleus-nucleus potential, the friction parameter, 
and the random force characterizing a one-dimensional Langevin-type description of the heavy-ion 
fusion process are investigated. Systematic calculations with the ImQMD model show that the 
fluctuation-dissipation relation found in the symmetric head-on fusion reactions at energies just 
above the Coulomb barrier fades out when the incident energy increases. It turns out that this 
dynamical change with increasing incident energy is caused by a specific behavior of the friction 
parameter which directly depends on the microscopic dynamical process, i.e., on how the collective 
energy of the relative motion is transferred into the intrinsic excitation energy. It is shown micro- 
scopically that the energy dissipation in the fusion process is governed by two mechanisms: One is 
caused by the nucleon exchanges between two fusing nuclei, and the other is due to a rearrangement 
of nucleons in the intrinsic system. The former mechanism monotonically increases the dissipative 
energy and shows a weak dependence on the incident energy, while the latter depends on both the 
relative distance between two fusing nuclei and the incident energy. It is shown that the latter 
mechanism is responsible for the energy dependence of the fusion potential and explains the fading 
out of the fluctuation-dissipation relation. 


PACS numbers: 24.60.-k, 24.10.Lx, 25.60.Pj, 25.70.Lm 


I. INTRODUCTION 


Energy dissipation process observed in heavy-ion reac- 
tions ranging from deep-inelastic collisions to sub-barrier 
fusions has been described by various macroscopic trans- 
port models [1-14]. Although many of these macroscopic 
models are very successful in predicting the cross section 
of compound nucleus formation, their microscopic foun- 
dation, e.g., how the colliding two nuclei fuse and how 
the relative kinetic energy dissipates into the intrinsic 
energy, still requires further research. 

Besides these macroscopic transport models, many mi- 
croscopic approaches like the time-dependent Hartree- 
Fock (TDHF) theory [15-28], the many-body correlation 
transport theory [29], and the quantum molecular dy- 
namics (QMD) [30, 31], the antisymmetrized molecular 
dynamics [32, 33], and the fermion molecular dynam- 
ics [34] models have been developed. It is worthwhile 
mentioning that, among the TDHF approaches, there 
have been proposed new methods in recent years like 
the dissipative-dynamics (DD) TDHF [22, 23] and the 
density-constrained (DC) TDHF [35, 36] theories, which 
intend to explore how the macroscopic collective behav- 
ior of two colliding nuclei described by the Langevin-type 
equation appears as a result of huge dimensional micro- 
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scopic dynamics. With the aid of these methods, the in- 
cident energy dependences of the macroscopic potential, 
the collective inertia parameter as well as the friction 
force have been derived from microscopic TDHF simula- 
tions and the discussions of which are made in relation 
with the internal structure change of the colliding nuclei. 

As is widely accepted, the above subject of heavy-ion 
fusion process is related to a more general problem, i.e., 
how to get a better understanding on various macro- 
scopic collective behavior of non-equilibrium systems in 
a deeper level, which is currently studied in many fields 
of sciences, e.g., in physical and chemical as well as bi- 
ological systems, by exploiting state-of-the-art ab-initio 
numerical simulations of the molecular dynamics model 
[37, 38]. 

Recently, we have proposed a macroscopic reduction 
procedure based on the improved quantum molecular dy- 
namics (ImQMD) model, aiming at exploring how the 
macroscopic Langevin-type equation emerges out of the 
microscopic dynamics in such a finite quantum many- 
body system as an atomic nucleus [39]. We found that 
the dissipation dynamics of the relative motion between 
two fusing nuclei is caused by a non-Gaussian distribu- 
tion of the random force. A clear non-Markovian effect 
was also observed in the time correlation function of the 
random force. Note that the non-Gaussian fluctuation 
is a general feature in non-equilibrium systems ranging 
from the cosmoscopic [40, 41], mesoscopic [42, 43], to 
microscopic systems [44]. Furthermore, as discussed in 
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Ref. [45], the non-Markovian dynamics were relevant in 
many fields and applications. 

In this paper, we further develop our macroscopic re- 
duction procedure applicable to the ImQMD model. An 
extension is based on a new method [46] which projects 
the effects of intrinsic degrees of freedom onto the col- 
lective subspace by transforming the time-scale into a 
macroscopic collective space-scale. With the aid of this 
macroscopic scale transformation, we explicitly extract 
the effects of intrinsic degrees of freedom in a convenient 
form that is appropriate for discussing the dynamical 
change of the fusion reaction from the adiabatic regime 
to the diabatic one. Exploiting the representation ob- 
tained after the transformation, one may clearly discuss 
how the dynamical role of the intrinsic system changes 
when the incident energy increases and how the energy 
transfer from collective motion to the intrinsic one, which 
is found to be carried out through both nucleon transfer 
between two fusing nuclei and the rearrangement effects 
in the intrinsic system. 

The paper is organized as follows. In Sec. II, we briefly 
introduce the ImQMD model and recapitulate the macro- 
scopic reduction procedure of the ImQMD model pro- 
posed in Ref. [39]. Applying the macroscopic reduction 
method of the ImQMD model to °°Zr+°°Zr head-on fu- 
sion reaction, in Sec. II, numerical results on the incident 
energy dependences of the random force, the fluctuation- 
dissipation relation, and the fusion potential are dis- 
cussed. In Sec. IV, we introduce various macroscopic 
quantities by further developing the macroscopic reduc- 
tion procedure suitable for exploring the macroscopic dis- 
sipative motion. Applying these new quantities on the 
numerical simulations of the InQMD model, it is clearly 
shown that the energy dissipation is characterized by two 
competitive microscopic processes. The energy depen- 
dence of the fluctuation-dissipation relation and that of 
the fusion potential are consistently explained as a re- 
sult of these competitive factors. Concluding remarks are 
given in Sec. V. In the appendix, we give the derivation 
of the energy dissipation due to the nucleon exchange. 


Il. IMQMD MODEL AND THE MACROSCOPIC 
REDUCTION PROCEDURE 


A. ImQMD 


Like in the original QMD model [30] and in various 
modern versions of QMD models [47-54], a trial wave 
function for a nucleus in the ImQMD model [55-59] is 
restricted within a parameter space {rj, pj}, where rj 
and p; are mean values of position and momentum op- 
erators of the jth nucleon expressed by a Gaussian wave 
packet, and the total wave function is a direct product 
of these wave functions of Gaussian form. With the aid 
of a Wigner transformation, a nucleus composed of dis- 
tinguishable N nucleons is characterized by the following 


one-body phase space distribution function, 
1 (r—r;)? 20? P 
fee) = gpa E- Sepp.) |) 


202 h 
and a density distribution function is expressed as 


where g, represents the spatial width of the nucleon wave 
packet. 

The time evolution of the system in question is gov- 
erned by a set of canonical equations of motion 


oF = (3) 


which is derived from the time-dependent variational 
principle [30, 32, 34]. Upon that, collisions are included 
in the ImQMD simulations and the scattering angle of 
a single nucleon-nucleon collision is randomly chosen in 
such a way that the distribution of the scattering angles 
of all collisions agrees with the measured angular distri- 
bution for elastic and inelastic collisions [30]. The Hamil- 
tonian H consists of the kinetic energy and an effective 
interaction potential energy, 


AH =T + Uc + Ucou- (4) 


Here, Ucou denotes the Coulomb energy and the nuclear 
potential energy Ujoc is expressed as 


Uoc Ea 
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In the present paper, Vioc(p| is defined by applying the 
effective Skyrme interaction energy density functional 
without the spin-orbit term. The isospin asymmetry 
ô = (Pn — Pp)/(Pn + Pp) where p, Pn, and pp are the 
nucleon, neutron, and proton densities, respectively. In 
the present work, the parameter set IQ2 is used [60]. 

After integration, one obtains the local interaction po- 
tential energy: 
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where 
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and ti = +1 for protons and —1 for neutrons. The 


Coulomb energy is expressed as the sum of the direct 
and the a contribution 


Ucoul = 


For low energy nuclear reactions, the Pauli princi- 
ple plays an important role [34, 61]. In the ImQMD 
model, by using the phase space occupation constraint 
method [62], the fermionic properties of nucleons is ap- 
proximately taken into account. 

Many ImQMD simulations are made and each of them 
is called one event. We examine the average properties of 
these events and the deviation of each event from the av- 
erage. For each event, to prepare the initial nuclei for the 
projectile and the target, the position and momentum of 
each nucleon are randomly given under certain macro- 
scopic conditions, such as the binding energy and the 
radius. Numerical details can be found in Refs. [63, 64]. 


B. Macroscopic Reduction Procedure for InQMD 
model 


Since we focus on symmetric fusion reactions with 
the impact parameter equal to zero, we can introduce 
a separation plane at the center of mass (CoM) of the 
whole system, and divide the whole system into the left- 
(projectile-like) and right- (target-like) half parts. The 
relative motion between two CoMs of the left and right 
parts is chosen to be the relevant degree of freedom to be 
described by the Langevin equation. 

The one-dimensional generalized Langevin equation 
with memory effects is given as [2, 65-67] 


dP(t) ' / f / 
-f ye- npea ara- FP) 
where P(t) is the relative momentum between two parts, 
OF (t) is the random force felt by either part, and U (R) is 
the collective potential between two parts. The first term 
on the right hand side represents the retarded friction 
force. 

Based on this stochastic equation, in the ImQMD sim- 
ulations, the mean value of the relative momentum (P) pr 
between two CoMs at a given R is defined as 


= FD PO) 


Vetsei(t=R}? (10) 


where P;(t) and R;(t) are the momentum and position of 
the i-th event at time t. The time t; for the i-th event 
is chosen in such a way that the relative distance takes 
a given value R, i.e., R;(t;) = R. In the present work, 
various microscopic variables are discussed as a function 
of relative distance R instead of time t, because R char- 
acterizes how near the two nuclei locate. In Eq. (10) and 
hereafter, (Q) denotes an average of Q over all events. 
When it does not cause any confusion, we hereafter use 
the same notation Q(R) for the single event Q;(R) as 
well as for the average (Q(R)) = S07_, Qi(R)/n other- 
wise mentioned. 

A collective potential for the relative motion is defined 
as, 


U(R) = Erot (R) — Nett (R) m Fright (R), (11) 


where Frot(R), Fier (R), and Erignt (R) represent the en- 
ergy of the total system and those of the left and right 
parts, respectively. Each of them consists of the kinetic 
energy, the nuclear and the Coulomb potential energies, 
and numerical results of U (R) are shown in Figs. 2 and 
3(b). The DC- and DD-TDHF have been applied to ex- 
tract the collective potentials between two nuclei [22, 68] 
which gave similar features as those obtained in the 
present work and in other QMD simulations [69, 70]. 

With the collective potential and the momentum, we 
define the collective energy as 


Econ (R) = Teou(R) + U(R), 


2 
Teou(R) = ile, (12) 
2u 
where u is the reduced mass of the system. After the 
two nuclei contact each other, u becomes dependent on 
the relative distance R and Eec.m. [22, 23, 64]. As our 
discussion is limited to regions of R and Ee.m. where pu 
does not change largely, in this work we take it to be a 
constant. 
In the ImQMD simulations, the random force for the 
i-th event is defined as 


oF; (R) 


A 
PRD RO 


= F'(R) - (F(R), (13) 


(14) 


{t|R(t)=R} ” 
: > F(R) (15) 
m i=1 


where f? (t) denotes, in the i-th event, the force acting 
on the j-th nucleon due to all other nucleons in the two 
fusing nuclei, A is the number of nucleons contained in 
the left (right) part, and n is the total number of events. 

The fluctuation-dissipation relation links the energy 
dissipation of the collective energy E.ou(R) to the fluc- 
tuation force originated from the microscopic dynamics, 
both of which play decisive roles in the macroscopic de- 
scription of dissipation phenomena. In the ImnQMD sim- 
ulations, the strength of the fluctuation is characterized 


= 
3 
Il 


201707.00921v1 


chinaXiv 


by 
(5F(R)SF(R)) = L S5 5F(R)SP:(R). (16) 


Assuming the work done by the collective motion 
against the friction force is completely converted into the 
intrinsic energy Fintr(R), we get a relation 

Eintr(R) = en, E coll (R), (17) 
where E..m, means the initial bombarding energy. With 
the aid of the Rayleigh’s dissipation function [71, 72] de- 
fined as 


F = p(R)(dR/dt}’, (18) 


where q( R) expresses the friction parameter without con- 
sidering the non-Markovian effects, the rate of energy 
loss from the collective motion (which is equivalent to 
the energy gain of the intrinsic system under the above 
assumption) is expressed as dEintr = Fdt. From these 
relations, one gets 


Ea) = (RS = 7(R)P. (19) 


The friction parameter yo under the Markovian approx- 
imation is expressed as 


= (Fi fric (R)) 

(R) = (PR , 

It should be noticed that the R-dependence of the friction 

parameter y(R) contains a sensitive information on the 

dynamics of energy transfer from the collective motion 

to the intrinsic degrees of freedom, because it depends 

on the derivative of the collective potential U(R) with 
respect to R. 


Frric(R) =——— m (20) 


II. INCIDENT ENERGY DEPENDENCE OF 
FLUCTUATION AND DISSIPATION 


A. Random force in the ImQMD simulations 


Figure 1 shows the width of the random force, which 
is defined as the root mean square of F(R), i.e., the 
square root value of the strength of fluctuation defined 
in Eq. (16). In the case of Eoin, = 195 MeV, the width 
of the random force turns out to be a smooth function of 
relative distance R, reaching a peak at R = 11 fm which 
locates slightly inside of the Coulomb barrier and then 
levels out after a small down slope. At R = 30 fm where 
we start our numerical simulation, the fluctuation comes 
from the randomness of the position and momentum of 
each particle at the initial time when each event is initial- 
ized under the macroscopic conditions. This fluctuation 
propagates following the equation of motion (3). Since 
the scattering angles of the nucleon-nucleon collisions are 


randomly chosen in such a way that the distribution of 
the scattering angles of all collisions agree with the mea- 
sured angular distribution for elastic and inelastic colli- 
sions in the QMD model [30], the randomness coming 
from the scattering may also contribute to the fluctua- 
tion. According to the shape of the distribution of the 
random force, we divide the whole process into three re- 
gions [39]. As is seen from Fig. 1, the width is narrow 
and stays unchanged at 15 fm < R < 30 fm (region J), 
indicating the stability of ImQMD simulations. 

In region II (11 fm Ș R < 15 fm) where the two nuclei 
are going to merge into one nucleus, one may see that the 
width has a growing asymmetric component in addition 
to the symmetric Gaussian. It turns out [39] that the 
asymmetric distribution is caused by a small number of 
events where a few nucleons are exchanged between two 
nuclei; these nucleons play a role in opening a window. In 
the majority of events which contribute to the main part 
of the Gaussian distribution, all nucleons well split into 
two separated groups, expressing a projectile nucleus and 
a target nucleus and keeping their initial stable mean- 
fields. As can be seen in Fig. 1, after eliminating those 
events where a few nucleons are exchanged between two 
nuclei, one gets a Gaussian distribution for the random 
force (the dark blue part in the inset corresponding to 
R = 13.5 fm) and a smaller root mean square of the 
random force displayed by the pink diamond in region 
Il. 

After this merging stage, we may see that the width of 
the random force in region ITI has again a Gaussian shape 
which is, however, two order of magnitude wider than 
that in region I. The main origin of this enlargement is 
caused by the above discussed transferred nucleons that 
feel the nuclear force from the other nucleus. It is an open 
question whether this factor is related to the formation 
of a neck or not [4, 8, 9, 12, 73-79]. 

Figure 2 shows how the strength of random force and 
the friction parameter yo(R) defined in Eq. (20) change 
depending on the initial bombarding energy Ee.m.. One 
may see that the fluctuation (F (R)dF(R)) and the fric- 
tion coefficient yo(R) at Eem. = 195 MeV show sim- 
ilar shapes and their peaks locate at the same point. 
Such a bumped shape in the friction parameter near 
the Coulomb barrier energy is also observed in the DD- 
TDHF calculations [23, 24]. 

When Eec.m. increases, the friction parameter exhibits 
a sizable energy dependence. Comparing the case of 
Ecm. = 195 MeV with that of Eem. = 205 MeV 
[Figs. 2(a) and (b)], one may recognize that the peak 
in the friction parameter disappears, while the strength 
of the random force keeps its shape almost unchanged. 
That is, a similarity between these two curves at ener- 
gies just above the Coulomb barrier gradually disappears 
as E..m. increases. When Ee.m. increases more, the fric- 
tion parameter shows a monotonically increasing behav- 
ior [Figs. 2(c) and (d)]. At the region near the barrier 
top 11 fms R < 12 fm, that is, the initial stage of fusion 
just after the capture, the higher the incident energy, the 
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FIG. 1. (Color online) Root mean square of the random force distribution as a function of relative distance R. Each inset 
shows a typical distribution in regions I, II, and III respectively. In the inset corresponding to a typical distribution in region 
II, the distribution is shifted so that the symmetric Gaussian distribution centers at ôF = 0, but the mean value of the random 
force is still zero. The pink diamond is the root mean square of the random force after eliminating events in the asymmetric 
tail as shown by the inset in region II. 
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FIG. 2. (Color online) The strength of random force (5F(R)dF(R))'* [red dots, in (MeV/fm)?] and the friction coefficient 
yo(R) (blue squares, in 0.001 c/fm) at Eom. = 195 MeV, 205 MeV, 215 MeV, and 235 MeV, respectively. The grey line shows 
the potential U(R). Pink diamonds represent (6F'(R)dF(R))**™ which is obtained by eliminating events in the asymmetric tail. 
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smaller the friction parameter, while the situation be- 
comes reverse when the two colliding nuclei come closer 
with each other. 

Note that in the ImQMD simulations, the nucleon col- 
lision effects are included in addition to the one-body dis- 
sipation. From the above discussions, it is seen that the 
ImQMD model gives similar energy dependence and mag- 
nitude of the friction parameter as those of the TDHF 
calculations [23, 24, 80], indicating a dominance of one- 
body dissipation in the energy region under discussion. 


B. Incident Energy Dependences of Collective 
Kinetic, Potential, and Intrinsic Energies 


Recently, it becomes clear that the collective potential 
U(R) exhibits rather sensitive energy dependence when 
Ec.m. takes near the Coulomb barrier energy [22, 80, 81]. 
It has been discussed that this energy dependence might 
be related to the dynamical change from adiabatic to sud- 
den fusion process: The two fusing nuclei have enough 
time to rearrange their intrinsic structure when Ee.m. 
is just above the Coulomb barrier energy, and the non- 
adiabatic effects gradually play a role as Ee.m. increases. 
However, the above arguments are mainly based on the 
bulk information like the density distribution of many 
nucleons, without referring to any microscopic dynam- 
ics of individual nucleons. Next we investigate the inci- 
dent energy dependences of the collective kinetic energy 
Teou(R) = Econ(R) — U(R), the potential energy U(R), 
and the intrinsic energy Eintr(R) with the microscopic 
ImQMD model. 

Figure 3(a) shows the R-dependence of collective ki- 
netic energy T-on(R) as a function of relative distance at 
Eom, = 195 MeV and Eom. = 235 MeV. In the case of 
Eom. = 195 MeV, the main part of the collective kinetic 
energy has already changed into a form of collective po- 
tential energy just after the collective motion overcomes 
the Coulomb barrier at R = 12.5 fm. Since there re- 
mains about 3 MeV rather constant collective kinetic en- 
ergy in 8 fm < R < 11.5 fm as is shown in the inset 
of Fig. 3(a), the fusing process takes place very slowly. 
In the case of Eom, = 235 MeV, after passing through 
the barrier, the collective kinetic energy takes almost ten 
times larger value than that in the case of Eom, = 195 
MeV, and shows continuous down slope without reach- 
ing a constant value. Here, it should be mentioned that 
the dissipated energy (equivalent to Eintr) at Eom, = 235 
MeV is much smaller than that at Eom, = 195 MeV in 
the region 9 fm < R < 11.5 fm as shown in Fig. 3(c), 
which is caused by the same reason why the friction pa- 
rameter gets reduced as Eem. increases. This point will 
be clarified in Sec. IV. 

By comparing Eintr(R) in Fig. 3(c) with U(R) and 
Teou(R) in Figs. 3(a) and 3(b), one finds that the dis- 
sipated energy at Ee.m. = 235 MeV is around 30 MeV 
smaller than that at 195 MeV, while U(R) is around 30 
MeV larger at R = 10 fm. This means that the friction 
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FIG. 3. (Color online) (a) Collective kinetic energy Teon( R) = 
Econ( R) — U (R), (b) potential energy U (R), and (c) intrinsic 
energy Eintr(R) as functions of relative distance R. The red 
line represents results at Ece.m. = 195 MeV and the blue line 
shows results at Hom. = 195 MeV. In (a), the inset shows the 
kinetic energy from 8 fm to 10 fm at Eom. = 195 MeV. In 
(b), the light grey line shows the potential energy calculated 
under the frozen approximation. 


force is not proportional but inversely proportional to the 
momentum of the collective motion at the initial stage of 
fusion process at R S 10 fm. In Fig. 3(b), the grey line 
shows the potential energy calculated under the frozen 
approximation, i.e., the potential calculated by fixing the 
density distribution of the projectile and target and mak- 
ing them overlapped at a given R. The barrier height is 
about 182 MeV. The deviation of U(R) from this po- 
tential energy under the frozen approximation could be 
attributed to the rearrangement effect of all the intrinsic 
degrees of freedom, which will be discussed in the next 
section. 
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IV. MICROSCOPIC DYNAMICS OF ENERGY 
DISSIPATION 


In the previous Section, we have discussed Ee.m.- 
dependence of the nucleus-nucleus potential U (R) as well 
as that of the friction parameter yo(J2). Since the signif- 
icant E..m.-dependence occurs near the Coulomb barrier 
energy where the applicability of the adiabatic approxi- 
mation is gradually replaced by that of the sudden ap- 
proximation when FE.m. increases, let us introduce new 
macroscopic quantities and further develop the macro- 
scopic reduction procedure. 


A. New Quantities Based on Macroscopic 
Reduction Procedure 


In this Subsection, we discuss how the collective en- 
ergy is transferred into the intrinsic one by dividing the 
dissipation into two processes. The first one is caused 
by the nucleon transfer which takes place frequently be- 
tween the left and right parts of the whole system, and 
the amount of energy transferred from the relative kinetic 
energy to the intrinsic one up to a point R is given by 
Taiss(R) in Eq. (A2) in Appendix A. 

To discuss the second process of dissipation, we pay 
attention to Fj,,(R) for the i-th event which expresses 
the difference between the force obtained by differentiat- 
ing the nucleus-nucleus potential U(R) in Eq. (11) and 
that obtained by summing all the nucleon-nucleon forces 
in Eq. (14). Ff e(R) is defined as, 


Fits (R) = Fis (R) — noun) 


F(R) =-" 1) rR) = F(R), 21) 


where U*(R) and F"(R) for the i-th event are given in 
Eqs. (11) and (14), respectively. For the sake of con- 
venience, Fi (R) and Fi,(R) will be called as the col- 
lective force and the total nucleus-nucleus force, respec- 
tively. Corresponding to the collective force F(R), one 
may introduce another potential U*.,,(R) defined by 


coll 


R 
Uin(R) = — f i n(R)dR', Uky(co) = 0, (22) 


which expresses the work done by the collective system 
against F(R) up to the point R. Using Ucon(R) to- 
gether with U (R) in Eq. (11), in addition to the differ- 
ence in Eq. (21), one may further introduce a deference 
as 


R 
W(R)=- f dR! (Fi, — Fiy) 


R 


oO 


As is well known, the reduction of many-body dynam- 
ics onto a one-dimensional collective space inevitably in- 
troduces additional effects coming from the intrinsic de- 
grees of freedom, besides the genuine effect within the col- 
lective space. That is, F(R) and W (R) originally act- 
ing in the intrinsic space are expressed as the additional 
force and additional potential in the collective space be- 
cause of eliminating the intrinsic degrees of freedom, and 
of changing the time-scale into the collective space-scale 
introduced in Eq. (10). More generally, it is shown an- 
alytically [46] that the dynamics in the collective system 
and that in the intrinsic system are distinguishable when 
one applies the macroscopic reduction procedure on the 
ImQMD simulations. 

The additional potential W(R) is then regarded as the 
work done by the intrinsic system to rearrange its state 
which has been disturbed by the transferred nucleons. 
In other words, —W (R) represents the second process 
which changes the energy in the intrinsic system. With 
the aid of W (R) and Taiss( R), the total amount of energy 
transferred from the relative motion to the intrinsic one 
is given by 


Eäiss (R) = Taiss (R) E W(R), (24) 


where the former is from the collective kinetic energy 
through the nucleon transfer, and the latter is due to the 
rearrangement of the intrinsic system. 

To test our calculation, we may pay attention to a new 
quantity defined through 

2 
= t UeanlR)}. (25) 
2u 


where Ucon( R) is defined in Eq. (22). It should be noticed 
that there holds an energy conservation law within the 
collective degree of freedom as 


Kaiss(R) = Eom. { 


2 


P 
Eom. = Taiss(R) + p 


Ju + Uon(R) } F (26) 


where Taiss( R) means the dissipated energy due to nu- 
cleon exchange defined in Eq. (A2) as mentioned at the 
beginning of this Subsection, provided that there holds 
the relation 


Taiss (2) = Kaiss(R). (27) 


In the following, the validity of Eq. (27) will be shown 
with numerical simulations. By using Eq. (24), the ap- 
proximate energy conservation law in Eq. (26) is then 
expressed as 


2 


P 
Eom = Edis a 
m. d Cuy 


lr Ucou(R) + W(R). (28) 
Since the relation 
Eintr(R) + Eaiss(R), (29) 


is satisfied [c.f. Fig. 5 and relevant discussions], Eq. (28) 
is reduced to Eq. (17) which expresses an energy conser- 
vation of the total system. 


(MeV/fm) 


intr 


F. 


FIG. 4. (Color online) Fintr(R), the difference between the 
total nucleus force Fot = —OU(R)/OR and the collective force 
Feou(R) for the case of (a) Ecm. = 195 MeV and (b) Eom. = 
235 MeV. “I” and “II” in each subfigure correspond to regions 
I and II defined in Fig. 1. 


B. Two Microscopic Processes of Energy 
Dissipation 


Making clear what information we get from the energy 
dependence of the collective potential U(R), we first ex- 
amine Fintr( R). It expresses what happens in the intrin- 
sic system when many nucleons are exchanged between 
two nuclei, and when two nuclei closely approach each 
other. It also tells us how these effects subsequently de- 
velop inside of the intrinsic system, which are usually 
considered as the rearrangement of the intrinsic system. 

As is seen from Fig. 4, our numerical simulation tells us 
that Fintr(R) takes positive value in the range of 10.5 fm< 
R< 12.5 fm (10.0 fm < R < 12.5 fm) and negative value 
at R < 10.5 fm (R <10 fm) in the case of Eem. = 195 
(Eom. = 235) MeV. It is also recognized from this figure 
that Fintr( R) in the case of Eom, = 195 MeV shows much 
stronger R-dependence than that of Fein, = 235 MeV. A 
steep slope of Fintr(R) in the region 9 fm < R < 11 fm 
at Eom, = 195 MeV indicates that the rearrangement of 
the intrinsic system takes place more strongly than that 
in the case of Eom. = 235 MeV. 

Since the intrinsic system is kept unchanged in the re- 
gion R > 12.0 fm irrespective of Ee.m., which is expected 
from Fig. 3(c), it is reasonable that Fintr(R) does not 
take any appreciable value in this region. In the case 
of Eom. = 195 MeV, the adiabatic approximation is ex- 
pected to be applicable because of a strong rearrange- 
ment of the intrinsic system. In the case of Eom, = 235 
MeV, however, there appears more gentle slope which 
indicates rather weak rearrangement effect in the intrin- 
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FIG. 5. (Color online) Two microscopic sources of dissipation 
and its comparison with an increase of intrinsic energy. (a) 
Taiss(R) is from nucleon exchange process. (b) W (R) is due to 
the rearrangement effects in the intrinsic system. (c) and (d) 
show the difference between two processes Faiss = Taiss( R) — 
W(R) compared to the intrinsic energy Fintr(R) at Hom. = 
195 MeV and Eec.m. = 235 MeV respectively. 


sic system. In this case, the main parts of original nu- 
clei invade into each other with much shorter distance R 
so that the frozen density and/or sudden approximation 
may have more sense than the case with Eem. = 195 
MeV. 

Let us discuss W(R) which is the integrated value of 
the difference Fintr(R). Equation (24) tells us that the 


increase of Egiss(R) comes from two sources, i.e., Tyiss(R) 
and W(R). As is discussed in Appendix A, the former 
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comes from the nucleon transfer between two fusing nu- 
clei, enhances the dissipation, and might be related to the 
window effect. On the other hand, the latter is caused 
by the rearrangement process of the intrinsic system. 
Namely, the disturbance due to the nucleon exchange, 
and the subsequent process of forming a new state like 
an excited mean-field in the intrinsic system play a role 
to first enhance the dissipation and subsequently to sup- 
press it. 

In Fig. 5, two microscopic sources of dissipation 
Taiss(R) and W(R), and their combined effects Egiss(R) 
in comparison with Eintr(R) defined in Eq. (17) are 
shown. Figures 5(c) and 5(d) show that there well holds 
the relation given in Eq. (29) irrespective of the initial 
bombarding energy. From this figure, one may conclude 
that the energy transfer from the collective motion to the 
intrinsic one in the region 8 fm < R < 11.5 fm can be un- 
derstood microscopically as follows. Firstly, the nucleon 
exchange process always contribute largely to the dissipa- 
tion and the dissipated energy Tuiss(R) can reach around 
150 MeV. Secondly, the rearrangement effect of intrinsic 
system may either enhance or suppress the dissipation, 
depending on the relative distance R. Thirdly, the rear- 
rangement effect is very sensitive to the incident energy. 
For the case of Eem = 195 MeV, the rearrangement ef- 
fect is very pronounced and —60 MeV < W(R) < 40 
MeV. This is because that, the fusion takes place slowly 
and there is much time for the whole system to rearrange 
nucleons. In other words, the fusion process goes adia- 
batically. However, for the case of Eem = 235 MeV, the 
rearrangement effect is less pronounced and —20 MeV 
< W(R) < 20 MeV. In this case, the fusion is mainly 
diabatic. The adiabatic or diabatic characteristics in the 
fusion will be discussed in more detail in next Subsection. 

From the above numerical results, one may see that the 
nucleus-nucleus interaction Ftot( R) is reduced to the sim- 
ple summation of constituent forces F.ou(2), provided 
there are no internal correlations in the intrinsic system. 
Namely Fintr(R) in Eq. (21) expresses the amount of cor- 
relations inside the intrinsic system. 

At the end of this Subsection, we pay attention to an- 
other approximate way of calculating Tyiss() given by 
Eq. (A2). In Fig. 6, we show the calculated Ucon(R) 
and Kaiss(R). Comparing Fig. 6(b) with Fig. 5(a), one 
can draw a conclusion that there holds the relation in 
Eq. (27). That is, the dissipation energy due to the nu- 
cleon exchange Tuiss(R) is well reproduced by Kgiss(R) 
defined in Eq. (25) which only takes account of the poten- 
tial Ucou(R) without considering the internal correlation 
W (R) of the intrinsic system. 


C. Adiabatic vs. Diabatic Processes 


In the previous Subsection, it is clearly shown that the 
fusion potential U (R) observed in Fig. 3(c) can be divided 
into two different components U.o1(R) and W(R) which 
are depicted in Figs. 6(a) and 5(b), respectively. Since 


in = 195 MeV 
m, 5235 MeV 
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FIG. 6. (Color online) (a) Calculated value of Ucou(R) as 
a function R. (b) Calculated value of Kaiss(R) = Eom. — 
P? /2u — Ucon(R) as a function R defined in Eq. (25). 


the energy dependence of W (R) is discussed in the pre- 
vious Subsection, let us explore the energy dependence 
of Taiss( R) which gives information on the incident en- 
ergy dependence of Ucon( R). From Figs. 5(a) and 6(b), 
the most important point to be explored on the energy 
dependence of Taiss(R) and Kaiss(R) is why the collec- 
tive energy dissipation through the nucleon transfer oc- 
curs more frequently at remote distance in the case of 
Eom, = 195 MeV than the case of Ee.m. = 235 MeV. 


In Fig. 7(a), the nucleon exchange rate defined as the 
number of nucleons which go through the separation 
plane per unit time is shown as a function of R. In the 
region 9.5 fm < R <13 fm, there are almost no apprecia- 
ble incident energy dependence in the nucleon exchange 
rate, though it becomes larger in the case of Eem. = 235 
MeV than the case of Eem. = 195 MeV for R < 9.5 
fm. In Fig. 7(b), the total number of nucleons Nx (R) 
which have ever transferred up to a distance R is shown. 
From Fig. 3(a), one may see that the kinetic energy of 
relative motion around R x 13 fm at Eom, = 235 MeV 
is more than 40 MeV larger than that at Eem. = 195 
MeV. Namely, the velocity in the former case is a few 
times faster than that in the case at Ecm. = 195 MeV. 
From Figs. 3(a) and 7(b), one may recognize that, at 
Eom. = 235 MeV, the nuclei approach each other too 
fast to exchange comparable amount of nucleons as in 
the case of Ee.m. = 195 MeV. In other words, the two 
nuclei at Eom, = 235 MeV keep much larger part of 
their original shapes into a small relative distance R than 
that at Eom, = 195 MeV. This mechanism explains the 
rather slow startup of the Tuiss(R) and Kaiss(R) in the 
case of Eem. = 235 MeV observed in Figs. 5(a) and 
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FIG. 7. (Color online) (a) Nucleon exchanging rate as a func- 
tion of relative distance, (b) Total number of nucleons which 
has ever transferred from left (or right) part to the right (or 
left). The red and blue lines represent two cases with Eom. = 
195 and 235 MeV, respectively. 


6(b). Since there holds the relation given in Eq. (25), 
the above discussed diabatic aspect of nucleon transfer 
process and the large collective momentum in the region 
8 fm S$ R < 11 fm for the case of Ee.m. = 235 MeV are 
two microscopic ingredients to determine the amount of 
work U.ou(R) done by the collective degrees of freedom. 


In order to show whether or not our microscopic discus- 
sion based on the role of Fintr( R) [and W(R)] is consistent 
with the understanding obtained by the density profile, 
the density at Eom, = 195 MeV and Ec.m. = 235 MeV, as 
well as that corresponding to the frozen density at differ- 
ent R are depicted in Fig. 8(a). The density profile under 
the frozen density approximation is obtained by fixing the 
density distribution of the projectile and target and mak- 
ing them overlapped at a given R by hands. An obvious 
difference in the density distribution among these three 
cases is observed at the central part or the “neck” region 
of the system. This situation is displayed more clearly by 
the densities along the reaction axis shown in Fig. 8(b), 
which shows the nucleon number densities along the reac- 
tion axis Z with X = Y = 0 at R= 10.5 fm. From these 
figures, one may learn that the density at Eom, = 195 
MeV in the neck region is obviously higher than the case 
at Eom. = 235 MeV, and the latter case gives a similar 
profile as that of the frozen density. This lower density 
at the neck region in the latter case just corresponds to 
the smaller number of transferred nucleons which explain 
the slower startup of the transfer process Taiss(R) in the 
case of Eom, = 235 MeV. 
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D. Friction parameter 


We are now in the position to discuss the energy de- 
pendence of the friction parameter yo(R) shown in Fig. 
2. According to the macroscopic reduction procedure, 
the friction parameter is derived in Eq. (20). Conse- 
quently, the energy dependence of the friction parameter 
is directly related to that of the nucleus-nucleus fusion 
potential which is shown in Fig. 3(b). From Figs. 3(b) 
and (c), one may recognize that the U(R) and Ejnty(R) 
at Eom. = 195 MeV show an inflection point at R = 
11.5 fm, while at Eom. = 235 MeV such an inflection 
point does not exist. This inflection point is responsible 
for the peak structure of the friction parameter shown in 
Fig. 2(a). 


Using Finty(R) ~ Eaiss(R), and from Eqs. (20) and 
(24), the friction parameter can also be expressed as 


vol R) = oe ( (T EY eo 


which tells us that the friction parameter is determined 
by two microscopic processes, i.e., the nucleon exchange 
and the rearrangement. Since the R dependence of 
Taiss(R) does not have the inflection point irrespective 
of the incident energy in Fig. 5(a), one may see that the 
first derivative of Tyiss(R) with respect to R does not 
show a strong incident energy dependence. On the other 
hand, —W(R) shown in Fig. 5(b) has a strong incident 
energy dependence. 


In Fig. 9, the R-dependence of the friction parameter 
is shown together with the two components originating 
from two different microscopic processes. The bumped 
structure of the friction parameter in the case of Eoin, = 
195 MeV can be understood to be a pronounced intrinsic 
effects at R ~ 11 fm [c-f. Fintr(2) shown in Fig. 4]. This 
bumped structure disappears rather quickly as Ee, in- 
creases, because the rearrangement effects of the intrinsic 
system around R = 11.0 fm become small, which is ob- 
served from Fig. 4. 


At the end of this Subsection, some comments on the 
validity of the fluctuation-dissipation relation should be 
addressed. As is mentioned above, the E..m, dependence 
of the friction parameter is directly related to the rear- 
rangement effects in the intrinsic system expressed by 
Fintr(R). On the other hand, the fluctuation 6F(R) is 
related to Feon( R) which does not depends on the rear- 
rangement taking place in the intrinsic system, and the 
shape change in the correlation function of the fluctua- 
tion does not show an energy dependence as is seen from 
Fig. 2. According to these different microscopic effects 
on the friction parameter and on the fluctuation force, the 
fluctuation-dissipation relation realized at Eom, = 195 
MeV breaks down when the incident energy increases. 
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FIG. 8. (Color online) (a) Density profiles obtained in ImQMD for different relative distances R at Fem.=195 MeV and 235 
MeV and those with the frozen density. The abscissa axis is the reaction axis Z and the vertical axis is axis Y. (b) Nucleon 
number density along the reaction axis Z with X = Y = 0 at relative R = 10.5 fm obtained by the ImQMD model. The red 
line represents the density at Ecm. = 195 MeV, the blue line at Fem. = 235 MeV. The dark grey dashed line is the density 
with frozen density, the light grey dashed-dotted lines are densities for projectile and target. 


V. CONCLUDING REMARKS 


In this paper, we have systematically studied the inci- 
dent energy dependence of the nucleus-nucleus potential, 
that of the friction parameter, as well as that of the ran- 
dom force in heavy-ion fusion reactions, by applying the 
macroscopic reduction procedure based on the ImQMD 
numerical simulation. Making the discussion clear, we 
focus our attention on the head-on fusion reaction of 
symmetric 9°Zr+°°Zr system at above the Coulomb bar- 
rier energies and pay special attention to the fluctuation- 
dissipation relation. 


It should be mentioned that our discussion is based on 


a separation of the total nucleus-nucleus force Fiot(R) 
into the collective part F.o1(R) and the intrinsic part 
Fintr(R). Although the latter also depends on the rela- 
tive distance R rather than the intrinsic coordinates, one 
may regard W(R) as the work performed by the intrin- 
sic system against Fintr(R). A theoretical justification 
of this point will be given in Ref. [46], where a reduc- 
tion of the many-body dynamics onto one-dimensional 
collective space is discussed. In the present paper, it is 
numerically shown that there holds the relation given in 
Eq. (27), i.e., Tuiss(R) ~ Kaiss(R) which gives another 
justification for the above separation. By exploiting the 
procedure appropriate for the ImnQMD numerical simu- 
lations, the first step toward the dynamics of emergence 
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FIG. 9. (Color online) R-dependence of friction parame- 
ter caused by two competitive microscopic processes at (a) 
Eom. = 195 MeV and (b) Fem. = 235 MeV. The blue dashed 
line shows effects of nucleon exchange. The red dashed line 
shows effects of rearrangement in intrinsic system. 


about how the macroscopic fusion dynamics comes out as 
a result of the microscopic nucleonic motion is clarified. 

We note that more study should be performed with 
non-zero impact parameters and for different reaction 
systems, particularly for asymmetric reactions. The spin- 
orbit coupling has been shown to be very important in 
low-energy heavy-ion fusion reactions [19, 27]. Further- 
more, it was also found that the negative shell correction 
energies lower potential barriers [53]. Therefore, more 
dissipations are expected if these effects are included in 
the ImQMD simulations. Finally, detailed investigations 
should also be made about the influence of two-body col- 
lisions on the dissipation process in fusion reactions at 
energies around the Coulomb barrier. 


Appendix A: Derivation of Taiss(R) 


Let us discuss the nucleon exchange process between 
two nuclei with mass M = Am where m denotes the nu- 
cleon mass and A the number of nucleons in the nucleus. 
Since we are interested in head-on symmetric heavy-ion 
collisions, we focus on the one-dimensional case depicted 
in Fig. 10 in which only the right part of the system is 
shown. 

At certain time to, the right part with mass mA is as- 
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FIG. 10. (Color online) The basic processes of nucleon ex- 
change between two fusing nuclei. 


sumed to move toward the left part with the momentum 
Po, and the kinetic energy T® = P? /2mA. In the fusion 
reaction, there occur many nucleon transfers between two 
nuclei which are characterized by two basic processes de- 
picted in Fig. 10. The one denoted as Process A is to 
get an additional nucleon with a momentum pı heading 
to the right, which is originally from the left part of the 
system. The other shown as Process B is to lose one nu- 
cleon with a momentum pı heading to the left. After the 
i-th nucleon transfer process, the kinetic energy of the 


right part denoted as TS ) is expressed as 


2 
TO — P; P; = Pi-1 — pi, Ai = Ai-1 + 7 
R WMA,’ i ? 


—1; nucleon from the right involved, 
, Ao = A, 


aa { +1; nucleon from the left involved, 


i=1,2,--- (Al) 
where 7 counts the number of nucleon transfer processes 
after to. The difference between the kinetic energy at the 
initial time to and that after the i-th nucleon transfer pro- 
cess is then expressed as Tuiss(i) = P?/2mA — P? /2mA;. 
We denote the number of nucleons which have been ex- 
changed up to a given relative distance R as Nex(R). 
Then one-half of the total dissipation energy of the col- 
lective motion at R is expressed as 


2 
PÈ PR) 


I iss R) = = 
i ( ) 2nA 2mAn.,.(R) 


ej 


Nex(R) 
Puan =Po- DS Po 
i=1 


Nex(R) 


Ana) = Aot X m (A2) 
i=1 


Here, it should be mentioned that the sum of the dissipa- 
tion energy in Eq. (A2) and the amount of increased en- 
ergy in the intrinsic kinetic energy of the right (left) part 
is easily proven to be zero. Namely, the amount of en- 
ergy lost from the collective motion through the nucleon 
exchange is just the same amount as that is increased in 
the intrinsic system. 
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